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SUMMARY 

A geometrically nonlinear finite-element analysis has been developed to 
calculate the strain energy released by delaminating plates during impact 
loading. Only the first mode of deformation, which is equivalent to static 
deflection, was treated. Both the impact loading and delamination in the 
plate were assumed to be axisymmetric. The strain energy release rate in 
peeling, Gj, and shear sliding, Gjj, modes were calculated using the frac- 
ture mechanics crack closure technique. Energy release rates for various 
delamination sizes and locations and for various plate configurations and 
materials were compared. The analysis indicated that shear sliding (Gjj) was 
the primary mode of delaraination growth. The analysis also indicated that the 
midplane (maximum transverse shear stress plane) delamination was more 
critical and would grow first before any other delamination of the same size 
near the midplane region. The delamination growth rate was higher (neutrally 
stable) for a low toughness (brittle) matrix and slower (stable) for high 
toughness matrix. The energy release rate in the peeling mode, Gj, for a 
near-surface delamination can be as high as 0.5 Gjj, and can contribute 
significantly to the delamination growth. 

INTRODUCTION 

Composite laminates made with high strength resin matrices are brittle 
and are easily damaged by transverse impact loads. Delaraination damage 
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commonly occurs under impacts by hard objects. Delaminations are usually 
invisible from the surface but still cause significant reduction in compres- 
sive residual strength of the laminate. The formation and growth of delamina- 
tions in laminated plates have recently been investigated under static loading 
(ref. 1). That study concluded that: (1) delaminations are formed in a 

laminate when the transverse shear around the contact area between the 
impactor and the plate exceeds a critical shear strength; (2) these delamina- 
tions grow if the shear strain energy release rate exceeds the matrix shear 
toughness. This report contains a finite-element formulation to predict the 
growth of an already created delaraination in a circular quasi-isotropic 
laminate using fracture mechanics criteria. The model delaminations shown in 
figure 1 were used to assess the effects of delamination size and location on 
the strain energy release rate and delamination growth. Although backface 
damage is one of the common failure modes in thin plates under impacts, that 
problem is not addressed here. Instead, possible delaminations in thick 
plates are examined. For simplicity, the impact event is assumed to be equiv- 
alent to an axisymmetric plate loaded slowly in the direction of impact. This 
assumption is reasonable for quasi-isotropic plates impacted at the center. 

Early studies at NASA Langley (refs. 2 and 3) of low velocity impact on 
8 to 32-ply circular laminates have indicated that laminates may undergo large 
deflections, which makes the problem geometrically nonlinear. The objectives 
of this study are: (1) to develop a geometrically nonlinear axisymmetric 

finite-element (F-E) analysis which includes the crack closure technique of 
calculating strain energy release rates, (2) to identify the mode (opening or 
shear sliding) of delamination growth due to low velocity impact, (3) to 
examine the influence of delamination size and location and plate config- 
uration on strain energy release rates, and (4) to examine the compression 


2 



Instability of the near-surface delaminated layer due to impact. Figure 1 
shows delamination models considered in this study. A midplane delamination 
(fig. 1(a)) was selected to identify the mode of delamination growth, i.e., 
mode l or 11 and/or combined. The critical location of the delamination was 
predicted by analyzing a plate with three delaminations; one at the midplane 
and the other two located at equal distances on each side of the midplane 
(see fig. 1(b)). 

Compressive flexural stresses that develop during the impact of a thick 
laminate may cause local instability of delaminated surface plies. The local 
buckling of the surface sublaminate results in high interlaminar peeling 
stresses in addition to shear stress at the delamination front. The mode of 
delamination growth therefore could be mode I, which is different from cases 
(a) and (b). To examine the energy release rates related to local insta- 
bility, a plate with a delamination nearer to the top surface of the laminate, 
referred to as a near-surface delamination, was analyzed. Although more than 
one delamination may be present or created during the impact, for simplicity 
they are neglected. 

A geometrically nonlinear two-dimensional finite-element program (ref. 4) 
was extended to include the axisymmetric case. The analysis includes calcula- 
tion of strain energy release rates in the two possible fracture modes, namely 
the peeling (G^) and the shear sliding (G-q) modes using the linear elastic 
fracture mechanics crack closure technique (ref. 5). Strain energy release 
rates for various delamination sizes and locations as well as plate configura- 
tions were calculated and compared. 
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SYMBOLS 


a plate radius, in 

F delamination front nodal force, N 

Gj strain energy release rate in peeling mode, J/nr 

2 

Gjj- strain energy release rate in shear sliding mode, J/m 

O 

G-£ c fracture toughness in peeling mode, J/nr 

G IIc fracture toughness in shear sliding mode, J/nr 

h plate thickness, m 

P static equivalent impact force, N 

q(r) contact pressure at r, N/m 4 

r-0-z plate cylindrical coordinate system 

R c impact pressure distribution (contact) radius, m 

Rg delamination radius, m 

u displacement in r-direction, m 

w displacement in z-direction, m 

w Q plate center deflection, m 

zg z-coordinate of a delaraination, m 

AR element size around the crack tip 

Subscripts: 

1,2 nodes behind the delamination-front node 

r,z radial and z-directions 

DESCRIPTION OF THE PROBLEM 

Figure 2 shows a circular, quasi-isotropic laminated plate of radius a 
and thickness h, loaded by an elastic sphere at the center. The cylindrical 
coordinates for the plate are r-9-z. A section at 9 = constant is shown in 
figure 2. The radial direction is r and the transverse direction is z. 

A midplane delamination of radius Rg is shown in figure 2. The z-coordinate 
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of the delamination is represented by Zgj for midplane delamination Zg =* 0. 
The impact by the sphere is assumed to be at the plate center and is repre- 
sented by a static equivalent force P. The deformation of the quasi- 
isotropic circular laminate and sublaminates in the delaminated region are 
assumed to be axisymmetric. The radial and transverse displacements of the 
plate are represented by u and w, respectively. The contact force between 

the sphere and the plate is represented by an elliptical pressure distribution 

f R c 

(ref. 6) (see fig. 3), such that P = 2 it I q(r)r dr, where q(r) is the 

•'0 

contact pressure at r. The radius R c of the pressure distribution depends 
on the impact force and stiffness and curvature properties of the sphere and 
the plate (see, for example, ref. 7). For the range of forces and the sphere 
and plate considered here, the value of R c varied from 0 to 1 mm. For 
simplicity, a nominal radius of 0.5 mm was assumed throughout the analysis. 

ANALYSIS 

A geometrically nonlinear, plane stress and plane strain finite-element 
(F-E) computer code GAMNAS (ref. 4) was extended to include axisymmetric 
analysis. The program uses a four noded isoparametric element. The method- 
ology followed in the analysis is given in reference 8. A finite-element 
model for a midplane delamination is shown in figure 3. The linear elastic 
fracture mechanics (LEFM) crack closure technique (ref. 5) was employed to 
calculate strain energy release rates in mode I and mode II at the delamina- 
tion front. A separate study was made to confirm the applicability of LEFM 
crack closure technique for calculating strain energy release rates in this 
geometrically nonlinear problem. Figure 4(b) shows a free- body diagram of the 
delamination front. Displacements in the r and z directions at nodes 1 
and 2 behind the delamination front (node 0) are, respectively, u-^, U£, w^. 
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and v *2 (see fig. 4(b)). Delamination tip nodal forces in the r and z 
directions are F_ and F z , respectively. Then, the strain energy release 


F Z I w 2 ” w ll 

rate in the peeling mode (I) is — ^T - 
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and in the shear sliding 


mode (II) is G 


F r l u 2 " u ll 

jj. = — jpjjg — — — , in which Rg is the delamination radius 
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and AR is the element size employed near the delamination tip. A fine F-E 
mesh (£R = 0.05 mm or 0.1 mm) was employed around the delamination front. 

Mesh refinement was studied to select a value of AR which yielded relative 
errors less than 1 percent. Strain energy release rates for various delamina- 
tion sizes, locations, and plate configurations were calculated and the 
results are presented in the following section. 


RESULTS AND DISCUSSION 

Strain energy release rates in mode I and mode II for the typical delami- 
nation configurations shown in figure 1 were calculated using the linear 
elastic fracture mechanics crack closure criteria. The plate was assumed to 
be loaded statically by an equivalent impact force acting over a radius R c 
(fig. 3). Energy release rates for various delamination sizes and location, 
and plate configurations and forces were calculated. Elastic material proper- 
ties used in the analysis are given in table I. Results are discussed in the 
following sections. 

A comparison of calculated Gj and G^-j for the delamination configura- 
tions shown in figure 1 indicated that for delaminations away from the sur- 
faces (figs. 1(a) and 1(b)) G-j- is very small, less than 0.1 percent of 
Gjj based on a deforraed-conf iguration analysis. Hence, Gj results are 
presented only for near surface delaminations. 
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Midplane Delamination 


Figure 5 shows the force P versus the plate center deflection w Q for 
eight delamination radii (0 to 12.7 mm). As expected, the increase in delami- 
nation radius reduces the plate flexural stiffness, hence the initial slope of 
the force-deflection curve decreases. 

, Figure 6 shows the variation of Gjj with force P for a range of 
delamination radii from 0.635 to 6.35 mm. An examination of the curves for 
various Rg shows that for a constant force P » G II decreases with an 
increase in delamination radius Rg > 1.27 mm. Therefore, the delamination 
growth stops after the initial growth. It is noted here that for R^ less 
than or equal to R c + h/2, the top compressive stresses reduce the effect of 
the shear stress singularity at the delamination front, with the result that 
calculated Gjj values were small (see dashed curve Rg = 0.635 mm). 

A method of predicting the delamination growth during loading is 
explained in figure 6. If a plate having mode II fracture toughness, say, 

Gj-£ C = 1.0 kJ/m^ and Rg = 1.27 (Rg > R c + h/2) is loaded from zero by a 
steel sphere, the Gjj increases from zero along the Rg = 1.27 curve 
(see fig. 6). At P = 0.87 kN, Gjj becomes equal to Gj-j c (1.0 kJ/m^). 
Further increase in the force P causes delamination to grow at constant 
Gji = Gjj c . The locus of Gjj is shown by a broken line in figure 6. The 
plot of Rg versus critical force P for a Gj Ic of a laminate is hereafter 
referred to as "delamination growth curve." 

Delamination growth curves for three materials whose elastic stiffness 
properties are the same as those for T300/5208 and with Gjj c values of 0.5, 
1.0, and 1.5 kJ/m^ are shown in figure 7. Each curve gives two values of 
Rg for a selected force P. The smaller value corresponds to the case where 
the delamination lies within the impact pressure region and the larger value 
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Is for outside the region. For < R c the compressive normal stress reduces 

the shear stress singularity at the de lamination front, hence larger force Is 
required for the delamlnatlon to become critical. This force could 

exceed the material strength In transverse shear, tension, and/or crushing 
(see ref. 3). Therefore, the creation and Initial growth of delamlnat Ions In 
and around the loaded region have to be governed by some other criteria like 
Interlaminar strength criteria. The horizontal dash-dot line shown In the 
figure Is the boundary (R^ > 1.1 mm) above which the present analysis (LEFM) 

Is valid. Failure In the region below the broken line Is controlled by 
Interlaminar strength criteria. A large delamlnatlon Rg > 1.27 ram grows 
raonotonlcally with the force. Furthermore, the rate of delamlnatlon growth 
(dR^/dP) is greater for low values of ^xic than f° r high values. Low values 
of ^xic designate low toughness material. 

The effect of impact force on energy release rate is shown in figure 8. 
Results for Rg = 4 mm and plate sizes 12.7 x 2 (radius x thickness in milli- 
meters), 12.7 x 3, and 25.4 x 2 are shown. Calculated decreased for 

both increased plate thickness (2 ram to 3 mm and radius 12.7 mm) and plate 
radius (12.7 mm to 25.4 mm) from the reference case (12.7 x 2). But, increas- 
ing thickness by 50 percent caused a much greater reduction in G-q than did 
doubling the plate's radius. 

Figure 9 shows G^^ versus force for T300/5208, C6000/HX205, C6000/F185, 
and T300/P1700 composites. Calculated G^'s f° r all four materials are 
practically the same, even though their stiffnesses differ (see table 1). 
However, a delamination in a brittle composite (like T300/5200), which will 
have lower Gjx c , will become critical and grow at lower force than in tougher 
composites. 
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Multiple Delaminations 

A laminate generally develops more than one delamination in thickness 
direction during the impact (see, for example, fig. 1(c)). A question is: 
which one of these delaminations is critical in the subsequent impact or 
transverse loading? To predict the location of a critical delamination, a 
plate with three delarainations was analyzed. One of the delaminations is at 
the raidplane (zg = 0) and among the remaining two, one is assumed to be above 
the midplane and the other below the midplane. Figure 10 shows strain energy 
release rate for three delaminations (Rg = 4 mm) located at Zg = 0 and 
±0.5 mm. The results indicate that Gjj for the midplane delamination is 
higher than that for the other two. Therefore, a delamination in the mid- 
plane, which is a plane of maximum transverse shear stress, is most critical 
and grows before other delaminations in that region. It may be noted that the 
present analysis is based on homogeneous laminate properties. In an actual 
laminate, maximum shear stress location depends on the stacking sequence which 
is a function of layup and polar angle 0 of the plate. 

Figure 11 shows the predicted critical forces for initiation of growth of 
three delaminations. The material fracture toughness Gjj c is assumed to be 
1.5 kJ/m . For small delaminations (Rg <2.5 mm) radii, the upper delamina- 
tion (A) grows first followed by B and C delaminations. However, for 
Rg > 3 mm, the midplane delamination grows first and followed by A and C. 
These curves can be used to predict delaraination radii for a given force by 
drawing an ordinate perpendicular to the abscissa as in figure 11. The inter- 
section of the ordinate and the curves gives the delaraination radii at the 
three locations. For example, an impact force of 2.5 kN causes delamination 
of radii 5.5, 4.4, and 4 mm at midplane, zg = -0.5, and 0.5 mm, respectively. 
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Near-Surface Delamination 


Compressive induced instability of delaminated laminas located near the 
impacted surface in a thick laminate was examined. If a thin debonded sub- 
laminate forms near the top surface of a laminated plate, it may buckle due to 
the compressive bending stresses developed by the transverse force on the 
plate (see fig. 1(c)). Large peel stresses caused by local buckling would 
tend to produce growth by mode I and mode II. Figure 12 shows energy release 
rates in both modes for delaminations located at zg = -0.75 mm at a force 
P = 2.4 kN. The two curves indicate that, due to local buckling, the ratio 
of Gj to Gjj could be as high as 0.5 (Rg = 4 mm). Since Gj c of a 
material is much smaller than G-j-j c , the delamination growth could be due to 
peeling alone. Both Gj and Gjj decrease with increase in Rg (Rg > 4 mm). 
Therefore, the delamination growth stops after an initial growth for a given 
force. For small values of Rg (Rg < 2 mm) Gj decreases and Gjj- 
increases. Although more than one delamination can exist in a laminated 
plate, for simplicity and to narrow attention to compression buckling, only a 
near surface delamination was considered in the above analysis. 

From figures 10 and 12 the values of G^j, for P = 2.4 kN and 
Rg * 4 mm, and for delaminations located at various values of Zg are given 
below. 


z g, mm 

-0.75 

-0.5 

0.0 

0.5 

Gjj, kJ/m 2 

0.21 

1.1 

1.4 

1.0 


These results indicate that the delaminations near the midplane would have 
reached critical G^j and grown before a near-surface delaraination becomes 
critical. The growth of near-surface delamination in mode II is, therefore, 
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less probable. However, the existence of mode, even though small, can 

cause surface delamination growth since a material's Gj c is far smaller than 

is G IIc* 

CONCLUSIONS 

An axisymmetric finite-element analysis of a quasi-isotropic circular 
laminate under a static equivalent impact force led to the following 
conclusions: 

1. The primary mode of delamination growth in the interior, away from 
top and bottom surfaces, of a circular, quasi-isotropic laminate subjected to 
transverse impact is by Gj-j (Gj-q was not addressed). 

2. For a constant force, G II decreases with increasing radius of the 
delamination when the delaraination radius is greater than the loaded region. 

3. Small delaminations (smaller than or nearly equal to the impact pres- 
sure area) have very low G-q's due to a reduced shear stress singularity. 
Therefore, initial interlaminar failures may be governed by the interlaminar 
strength criteria. 

4. Delamination growth rate is rapid (neutrally stable) for low fracture 
toughness (brittle) materials, and it is slow and stable for high toughness 
materials. 

5. A midplane delamination, with a radius larger than the impact pres- 
sure radius, is more critical and grows faster than other delaminations of the 
same size but located in any other plane. 

6. If a thin debonded sublarainate forms near the impact surface of a 
thick laminate, the sublaminate may buckle due to bending compressive 
stresses. The peel mode energy release rate G-j- can be as high as O.SGjj. 
Therefore, a near-surface delamination in a material with low G^ c may grow 
in the peeling mode. 
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Table 1. Elastic properties of materials 


Parameter 


Material 


T300/5208 

C6000/HX205 

C6000/F185 

T300/P1700 

E rr , GPa 

54.00 

48.33 

44.35 

37.56 

E zz’ GPa 

14.30 

9.17 

6.79 

6.64 

G rz , GPa 

3.50 

4.10 

2.70 

2.62 

v rz 

0.28 

0.27 

0.27 

0.38 

V r0 

0.31 

0.31 

0.32 

0.32 
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(a) Finite-element mesh around 
the delamination front. 


(b) Free-body diagram at the 
delamination front. 


Fig. 4 Nomenclature in the region of delamir.ation front 









Fig. 5 Impact force versus deflection for a quasi-isotropic, circular plate 
(a = 12.7 mm, h = 1.0 mm, T300/5208) " 
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Fig. 7 Critical force versus delamination radius for various values of 
(a - 12.7 mm, h = 1.0 mm, T300/5208) 
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Fig. 


Effect of force on for various materials for a mid-plane delamination plate 

(a = 12.7 mm, h = 1.0 mm, R. = 3.8 mm) 
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Fig. 11 Effect of force on radius of delamination for multiple delaminations, 
(a = 12.7 mm, h = 2.0 mm, Gj = 1.5 kj/m 2 , T300/5208) 
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Fig. 12 Strain energy release rates for near-surface delaminations for P 

(a = 12.7 mm, h = 2 mm, z = - 0.75 mm, T300/5208) 
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